Space-filling bearings in three dimensions 
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We present the first space-filling bearing in three dimensions. It is shown that a packing which 
contains only loops with even number of spheres can be constructed in a self-similar way and that it 
can act as a three dimensional bearing in which spheres can rotate without slip and with negligible 
torsion friction. 

PACS numbers: 46.55.-|-d, 45.70.-n, 61.43.Bn, 91.45.Dh 



Space-filling bearings have been introduced in several 
contexts, such as in explaining the so-called seismic gaps 
P, of geological faults in which two tectonic plates 
slide against each other with a friction much less than 
the expected one, without production of earthquakes or 
of any significant heat. Space- filling bearings have also 
been used as toy models for turbulence and can also be 
used in mechanical devices Two dimensional space- 
filling bearings have been shown to exist and a discrete 
infinity of realizations has been constructed [113. The 
remaining question still open is: Do they also exist in 
three dimensions? This question is of fundamental im- 
portance to the physical applications. 

In this letter, we will report the discovery of a self- 
similar space-filling bearing in which an arbitrary cho- 
sen sphere can rotate around any axis and all the other 
spheres rotate accordingly without any sliding and with 
negligible torsion friction. 

In two dimensions, different classes of space-fillinj 
bearings of disks have been constructed in Refs. ^ 
by requiring the loops to have an even number of disks, 
since in two dimensions this is obviously a necessary and 
sufficient condition for disks to be able to rotate without 
any slip. Successive disks must rotate, in alternation, 
clockwise and counter-clockwise in order to avoid frus- 
tration. 

The situation in three dimensions is different from two 
dimensions in two ways; The axes of rotation need not 
be parallel, and the centers of spheres in a loop may not 
lie all in the same plane. As a result, even in an isolated 
odd loop spheres could rotate without friction. But, as 
we will see, in the packings with an infinite number of 
interconnecting loops, we could construct unfrustrated 
configurations of rotating spheres when all loops have an 
even number of spheres. Such a packing is bichromatic, 
i.e., one can color the spheres using only two different 
colors in such a way that no spheres of same color touch 
each other, as shown in Fig^ 

No three dimensional space-filling bearing has been 
known up to now. The classical Apollonian packing is 
space-filling and self-similar but not a bearing since at 



FIG. 1: Bichromatic packing of spheres with octahedral sym- 
metry. No two spheres of same color touch each other. The 
image on the left shows the initial configuration and the first 
generation of inserted spheres. 



least five colors are needed to assure different colors at 
each contact. This packing can be constructed in differ- 
ent ways 0,0 ■ By generalizing the inversion technique 
used in Ref.ja] to other Platonic Solids than the tetrahe- 
dron (the base of 3D Apollonian packing) we were able 
to construct five new packings including a bichromatic 
one. Details on the construction algorithm and on the 
complete set of new configurations will be published else- 
where 0. We give here only a qualitative description of 
this technique for the bichromatic packing. 

Let us consider filling a sphere of unit radius. The 
filling procedure is initialized by placing seven initial 
spheres on the vertices and the center of a regular oc- 
tahedron inside the unit sphere, so that the spheres on 
the vertices do not touch each other but touch the unit 
sphere and the one in the center. Further spheres are 
inserted by inversion such that this topology can be 
preserved on all scales, imposing that, all spheres are on 
vertices or centers of (deformed) octahedra. 
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The inversions are made iteratively with respect to 
nine inversion svheres\l^: One inversion sphere is con- 
centric with the unit sphere, and is perpendicular to the 
six initial spheres on the vertices of the octahedron |Q . 
Inversion around it, therefore, leaves the vertex spheres 
invariant, and maps the unit sphere onto the central 
sphere. The other eight inversion spheres are associated 
to each face of the octahedron, that is, they are perpen- 
dicular to the unit sphere and to the three spheres on the 
vertices of that face. The inversion around each of these 
inversion spheres maps the unit sphere and spheres of 
the corresponding face into themselves (that is, it gives 
no new sphere), and the other four initial spheres into 
four new spheres within the space between the face and 
the unit sphere. 

Figure |2] shows the plane cut through the centers of 
the unit sphere and four initial spheres. Dashed circles 
are cuts of the inversion spheres. Sphere S : 0, 1, • • • is 
mapped, by the inversion sphere shown by a thick dashed 
line, onto sphere S'. The inversion around this sphere 
gives no new images of spheres 1 and 2. In the first it- 
eration we make all possible inversions which give new 
and smaller spheres. In the next iterations, the newly- 
generated spheres are mapped to smaller spheres. For 
example, sphere 0' is mapped (by the central inversion 
sphere) onto 0". In this way, the remaining empty space 
is filled in the limit of infinite iterations while the bichro- 
matic topology of the contacts is preserved. 

Using this algorithm, the configuration of initial 
spheres which gives a bichromatic packing is unique. 
Strictly speaking, it is shown that the only value for radii 
of the initial spheres which leads to the bichromatic pack- 
ing without (partial) overlapping of generated spheres is 
(•\/3 — l)/2 and only using an octahedral base. The frac- 
tal dimension of the packing has been computed using 
the same method as in Ref.jj] and it is 2.59 considerably 
higher than that of Apollonian packing in 3D, that is 
2.470. 

The image on the left of the FigQ] shows the initial 
configuration and the first generation of inserted spheres, 
and the one on the right shows the resulting packing con- 
taining all the spheres with radii greater than 2^^. The 
sphere at the center and the external hull are white and 
those on the vertices of the octahedron are black. Since 
the spheres on the vertices touch only the external hull 
and the central sphere, and since this topology is pre- 
served by construction, no two spheres of the same color 
touch each other. 

This implies that every loop of spheres in this packing 
contains an even number of spheres. We now show that 
this is a sufficient condition for the spheres in contact to 
rotate without slip, or even torsion friction. 

Consider a loop of n spheres as seen schematically for 
71 = 4 in FigI3 The no-slip condition implies that each 
pair of touching spheres have the same tangent velocities 
V at their contact point. The condition for the contact 




FIG. 2: Plane cut through the center of the unit sphere and 
four initial spheres. Dashed circles are cuts of the inversion 
spheres. Sphere S is mapped, by the inversion sphere shown 
by a thick dashed line, onto sphere 5". 0' is mapped by the 
central inversion sphere onto O". 

between the first and the second sphere can be written 
as: 

Vi = V2 

i?lfi2 X UJi = -R2fi2 X a;2 

{RlUJl + i?2'^2) X fi2 = 0, (1) 

where i?2, (^i and cj2 are the radii and the vectorial 
angular velocities of the first and second sphere, respec- 
tively. fi2 is the unit vector in the direction connect- 
ing the centers of the first and the second sphere. From 
Eq.(^) the vector {Riuii + i?2'^2) should be parallel to 

R2uj2 = -Riuji - ai2fi2, (2) 

where ai2 is an arbitrary parameter. Equation (j^J is a 
connection between the rotation vectors lui and uj2 of the 
two spheres in contact. Similarly for the third sphere in 
contact with the second, wc have 

= -i?2W2 - a23?^23- (3) 
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Putting Eq.|(21) into Eq.lO we find the relation between 
the angular velocities of the first and third sphere: 

R3UJ3, = RiUJi + ai2fi2 - a23?'23- (4) 

In general, we can relate the angular velocities of the first 
and jth spheres of an arbitrary chain of spheres in no-slip 
contacts by: 

i-i 

R,C3, = i-iy-^Riuji + ^(-iy-^a,,,+in,,+i. (5) 

i=l 

As long as the chain is open, the spheres can rotate with- 
out slip with the angular velocities given by Eq.© and 
no restrictions on a^.i+i. But, for a loop of n spheres in 
contact, spheres j and j + n are identical, so that 

n 

Ricji = + ^(-l)"-'+ia,,,+if,,,+i. (6) 

1=1 

A similar equation holds for every sphere j = 1, • • • , n in 
the loop. 

Although for a single loop there are many solutions 
of Eq.@, not all will serve our purpose. In a packing, 
each sphere belongs to a very large number of loops and 
all loops should be consistent and avoid frustration. In 
other words, the angular velocity obtained for a sphere 
as a member of one loop should be the same as being a 
member of any other loop. 

If the loop contains an even number n of spheres, 
Eq. © becomes a relation between the hitherto arbitrary 
coefficients of connection a^^i+i, 

n 

^(-l)^a,,,+if,,,+i = 0. (7) 
Using the fact that the loop is geometrically closed: 

n 

J2iR^+R^+l)f^.^+1^0, (8) 

a solution for Eq.© is 

a^,^+l=c{-ly{R^+R^+l), (9) 

where c is an arbitrary constant. Putting this in Eq.©, 
yields the angular velocities 

LS, = ^{~iy[-Ri^i + cR,,), (10) 

where Rij is the vector which connects the centers of 
the first and jth sphere. As can be seen, the angular 
velocities only depend on the positions of the spheres, 
so that the consistency between different loops can be 
automatically fulfilled providing that the parameter c is 
the same for every loop of the entire packing. 




FIG. 3: Schematic configuration of a loop of four spheres. 
The spheres rotate without slip. The centers of spheres need 
not be in the same plane, although cDi, uj2 and R12 are copla- 
nar. 



In Ea. (|10|l . all the angular velocities are calculated 
from LOi and c, which can be chosen arbitrarily, (c — 
corresponds to the case when all angular velocities are 
parallel.) 

The no-slip condition (0) then reads 

i?2W2 = C^12, (11) 

so that the vectors cDi, L02 and R12 are coplanar (the 
plane of Fig. 3, containing the two centers and the point 
of contact A). They are in general not coUinear. 

We note that the condition of rotation without slip 
(PI also guarantees that there will be no torsion friction, 
as long as the three vectors are not collinear. Indeed, 
the locus of the contact point A on sphere 1 is a circle 
Ci, perpendicular to cJi. On sphere 2, it is a circle C2, 
perpendicular to 002- The cone tangent to sphere i — 1,2 
on circle Ci has an apex Si, on the axis LOi. The two apices 
Si are in the same plane as the sphere centers and on the 
same line as the contact point A (Fig0J. The two spheres 
rolling on each other can therefore be replaced by the two 
cones rolling on each other, around the line S1S2A, which 
always contains the contact point A. There will be no 
twisting between cones, thus no torsion friction from one 
sphere rolling on the other. Only when the two angular 
velocities and R12 are collinear can there be some twist 
of the spheres against each other. (The circles Ci then 
reduce to the contact point A, there are no tangent cones, 
and the tangent velocities Vi = 0). This situation is not 
generic. It is of measure zero and physically irrelevant. 
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agrees with that we obtained in this work. 

In summary, we proved the existence of the three di- 
mensional space-fiUing sphere bearing by presenting the 
exphcit construction. We have shown that a sufficient 
condition is that the packing is bichromatic, and given 
an exphcit expression IjlOl) for the angular velocity of ev- 
ery sphere of the entire packing, in terms of wi and c 
only. In this way, we support the previous modclization 
for lubrication between tectonic plates. This result can 
also be important in mechanics and hydrodynamics. 

We acknowledge the hospitality of the Max Planck In- 
stitute for the Physics of Complex Systems, GEOMES 
program, oct. 2002. 



FIG. 4: Two spheres rolling on each other without slip can 
be represented by two cones rolling on each other around the 
common line S1S2A. Therefore there is no twist of the spheres 
against each other. 



In the case of odd loops, Eq.® becomes 

n 

^(-l)*a,,,+if, ,,+1 = 2Rii3i. (12) 

1=1 

A similar equation holds for every sphere j = 1, • • • , n in 
the loop. But, since the coefficients a^^i+i depend then on 
both geometry and the rotation velocities of the spheres 
of the loop, consistency between different loops cannot be 
fulfilled in general and, therefore, a packing containing 
odd loops cannot be a bearing. A rigorous proof for this, 
however, is missing. It should be mentioned however that 
an unfrustrated single odd loop is possible, but cannot 
occur in isolation in a packing [l(j| . 

The bearing discussed here is very idealized and based 
on exactly-spherical particles with infinite rigidity and, 
of course, does not exist in this form on all length scales 
in nature. Nevertheless, chances still remain that similar 
bearings, though with some imperfections, occur in real- 
ity. A simulation of two-dimensional shear bands shows 
formation of spontaneous rotating bearings in clusters of 
up to 30 particles Despite of having more volume, 
the bearing state is favored because of its low friction. As 
an another evidence, the bichromatic packing presented 
here is self-similar which is also observed in the sarnples 
of tectonic gouge down to several scales (see R,ef.|l2|V 
Interestingly, the measured fractal dimension, 2.60 ±0.1, 
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